APPENDIX A 



Determining Scanned Image First Axis Centerpoint of Reference 

int yellow_r=200 f yellow_g=200 / yellow_bl=50,yellow_b2=200; 

X=0; 

ReTryHoriz: 
X+=2 ; 

small c=bmp .bmHeight, bigc=0 ; 

for (y=bmp . bmHeight/8 ,*y<bmp .bmHeight-bmp .bmHeight/8 ;y++) 
{ 

GiveMeSome ( ) ; 
c=GetPixel (mDC,x,y) ; 

if(c!=white GetRValue <c) >yellow_r GetGValue (c) >yellow_g && GetBValue (c) >yellowJ 

GetBValue (c) <yellow_b2) 

{ 

if {y<smallc) 

smallc=y; 
if (y>bigc) 

bigc=y; 

} 

if (smallc==bmp. bmHeight && bigc==0 &S= x<100) 
{ 

yellow_r--5 ; 
yellow_g-=5; 
yellow_bl- =3 ; 
yellow_b2+=3 ; 
goto ReTryHoriz; 

} 

if (smallc==bmp.bmHeight && bigc==0) 
{ 

char str [100] ; 
SelectObject (mDC,old) ; 
DeleteObject (MyBmp) ; 
DeleteDC(mDC) ,- 
x=bmp . bmWidth ; 
y=bmp . bmHeight ; 

sprintf (str, "Can' t find horiz grid, please check prescan area <%lu, %lu) . " , x,y) ; 
MessageBox(hDlg,str,maintitle,MB_OK | M3_ICONSTOP) ; 

PostQuitMessage (FALSE) ; 

EndDialog (hDlg, FALSE) ; 

return (FALSE) ; 

} 

gridy=smallc+ (bigc-smallc) /2; 



Determining Scanned Image Second Axis Centerpoint of Reference 

int yellow_r=200,yellow_g=200 l yellow_bl=50,yellow_b2=200; 

y=0; 

ReTryVert : 
y+=2; 

small c=bmp .bmWidth, bigc=0 ; 

for <x=bmp .bmWidth/8 ; x<bmp . bmWidth-bmp . bmWidth/ S ; x++ ) 
{ 

GiveMeSome ( ) ; 
c=GetPixel{mDC,x,y) ; 

if(c!=white GetRValue (c) >yellow_r GetGValue (c) >yellow_g GetBValue (c) >yellow J 

GetBValue (c) <yellow_b2) 

{ 

if (x<smallc) 

smallc=x; 
if (x>bigc) 

bigc=x; 

} 

} 

if (smallc==bmp. bmWidth && bigc==0 && y<!00) 
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{ 

yellow_r-=5; 
yellow_g-=5; 
yellow_bl-=3 ; 
yellow_b2+=3 ; 
goto ReTryVert; 

} 

if (small c==bmp. bmWidth bigc==0) 

{ 

char str [200] ; 
SelectObject {mDC,old} ; 
DeleteObject (MyBmp) ; 

DeleteDC(mDC) ; 
x=bmp . bmWidth ,- 
y^bmp . bmHeight ; 

sprintf (str, "Can 1 1 find vertical grid, please check prescan area (%lu-%u, %u) 11 , gridy,x, y) ; 
MessageBox(hDlg,str,maintitle,MB_OK | MB_IC0NST0P) ; 

PostQuitMessage (FALSE) ; 

EndDialog (hDlg, FALSE) ; 

return (FALSE) ; 

} 

gridx=smallc+ (bigc-smallc) /2 ; 



Determining A Starting Radius 

BOOL bluecheck(COLORREF c) 

{ 

if ({unsigned char) Get RValue (c) < (unsigned char) 185 && (unsigned char) GetGValue (c) < (unsigned 
char) 185 && (unsigned char) GetBValue (c) > (unsigned char) 50) 
return (TRUE) ; 
return (FALSE) ; 

} 

if (bmp . bmHeight >bmp . bmWidth) 

{ 

if (bmp. bmHeight >3300) // 600dpi, 5.5 inches 

bmp .bmHeight =bmp .bmWidth ; 
g-bmp . bmHeight ; 

} 

else 

{ 

if <bmp.bmWidth>33 00) // 600dpi r 5.5 inches 

bmp . bmWidth=bmp . bmHeight ; 
g=bmp . bmWidth; 

} 

bigc2=bigc=0 ; 
smallc2=smallc=g; 

/////////// 

startmeof f =g; 

for (x=0;X<512;X+=128) 

{ 

double r; 

long int myx,myy; 

GiveMeSome { ) ; 

c=white; 

for (r~100 ; r<g && (c-=white || ibluecheck (c) ) ; r++) 

{ myx={int) ( (double) gridx+ (double) MySin (x) *r) ; 

myy=(int) ( (double) gridy+ (double) MyCos (x) *r) ; 

c=GetPixel (mDC,myx, myy) ; 

GiveMeSome ( ) ; 

} 

GiveMeSome () ; 
if {r<startmeof f ) 

startmeof f=r; 

} 

if (startmeof f >=g) 
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startmeof f =100 ; 

else 

{ 

startmeof f= (long) fabs (startmeof f/3) ; 
if (startmeof f <100) 

startmeof f =10 0 ; 

} 



Centering A Scanned Image Shape 

for (x= 0 ; x<MYPOINTS + 1 ; X++ ) 
{ 

double r; 

long int myx , myy ; 

double rads=x; 

c=white; 

for (r=startmeof f ;r<g (c==white || I bluecheck (c) ) ; r++} // our radius 

{ 

GiveMeSome ( ) ; 
if (reverseclick==0) 
rnyx=(int) ( (double) gridx- (double) MySin (rads) *r) ; 
else 

myx=(int) ( (double) gridx+ (double) MySin (rads) *r) ; 
tnyy= (int) ( (double) gridy+ (double) MyCos (rads) *r) ; 

c-GetPixel (mDCjmyx^yy) ; 

) 

if(r>=g) // means it's the max value.. 
{ 

char str (200] ; 
SelectObject (mDC,old) ; 
DeleteObject (MyBmp) ,* 
DeleteDC(mDC) ; 

sprintf (str, "No tracing found(%lu), please verify that your prsscan region is correct and that the 
tracing was done using an approved e.lens pen. n ,x); 
MessageBoxthDlg, str,maintitle,MB_OK | MB__IC0NST0P) ; 

Post Quit Message (FALSE) ; 

EndDialog (hDlg, FALSE) ; 

return (FALSE) ; 

} 

if (myx>bigc) 

bigc=myx; 
if (myx<smallc) 

smallc=myx; 
if (myy>bigc2) 

bigc2=myy; 
if (myy<smallc2) 

smallc2=myy; 
sprintf (instring, nu" , 40+30*x/MYPOINTS) ; 
SetDlgltemText (hDlg, BAR_CLICK, instring) ; 
InvalidateRgn(GetDlgItem(hDlg,BAR_CLICK) , NULL, TRUE) ; 
GiveMeSome () ; 

} 

gridx=smallc+ (bigc-smallc) /2; 
gridy=smallc2+ (bigc2-smallc2) /2 ; 



Determining A Scanned Image Radial Shape 

el ens . jobdata_datasize=0 ; 

for (x=0;X<MYPOINTS+l;X++) 

{ 

double r; 

long int myx, myy; 

double rads=x; 

c=white; 
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f or (r=startmeof f ;r<g (c==white || Ibluecheck (c) ) ; r++) 

{ 

GiveMeSotne () ,- 
if (reverseclick^=0) 
myx- (int) ( (double) gridx- (double) MySin (rads) *r) ; 
else 

myx= (int) ( (double) gridx+ (double)MySin(rads) *r) ; 
myy= (int) { (double) gridy+ (double) My Cos (rads) *r) ; 

c=GetPixel (mDC,myx,myy) ; 

} 

myx- (int)r; 
if (myx<l) 

myx-1 ; 

elens. jobdata_datasize++; 



Determining A Scanned Image Radial Size 
for (x*0;x<MYPOINTS+l;X++) { 

elens, shape [x] = ( ( (double) myx/ (double) config. resolution) / (double) 0 .0393 70) * (double) 100; 
elens * shape [x] = (unsigned short int) { (unsigned short int) ( (double) elens . shape [x] / (double) 10) * 10L 
+config. calibrate) ; 

} 



Smoothing A Scanned Image Radial Shape 

for (x*0;x<MYPOINTS+l;x++) { 

if {conf ig.makesmooth—1 && x>2) 

{ 

if (elens. shape Ix] >=elens. shape [x-2] ) 
{ 

if (elens. shape [x-1] <elens. shape [x-2] | | elens .shape [x-1] >elens. shape tx] ) 

elens. shape [x-l] = (short int) (elens. shape [x-2] + (elens. shape tx] -elens .shape [x-2] )/2) ; 
} 

} 

} 



Modify Size of Derived Radial Shape 

short int y; 
short int num=l; 
double val=CalcCircum() 
do 

{ 

for (y=0 ;y<5l3 ;y++) 
{ 

if (increasesize==l) 

elens , shape [y3 - {short int) (original [y] +of f set+num) ; 

else 

elens . shape [y] = (short int) {original [y] +of f set-num) ; 

} 

ShapeToJob ( ) ; 
JobTo Shape 0 ; 
num++ ; 

} while ( (increasesize==l &&val>=CalcCircum () ) || 
(increasesize==0 && val<=CalcCircum() ) ; 
if (increasesize==l) 

offsets (short int) (of f set+ (num-1) ) ,- 

else 

off set= (short int) (of f set- (num-1) ) ; 
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Rotating A Derived Radial and Smoothed Shape 

if (nasalup==l) 

{ 

short int a=original [0] ; 

for {y=0;y<elens. jobdata_datasize-l;y++) 

{ 

short int b; 
b=original [y+1] ; 
original [y+1] =a; 
a=b; 

} 

original [0] =a; 
rotated- - ; 

} 

else 

{ 

short int aboriginal [0] ; 

for (y=0 ;y<elens , jobdata_datasize-l ;y++) 

original [y] =original [y+1] ; 
original [y] =a; 
rotated++ ; 

} 
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